import scipy.special as sl
import numpy as np
import matplotlib.pyplot as plt

x = np.linspace(-10, 10, 1000)
y = np.linspace(-10, 10, 1000)
r = np.linspace(3, 5, 1000)
# r=np.sqrt(x**2+y**2)
z = sl.jv(0, r)
plt.figure()
plt.scatter(r, z, 0.1)

plt.show()

k = 1
while k <= 999:
    if z[k + 1] > z[k]:
        print(z, k)
        break
        pass
    k += 1
    pass
